Incomplete reprogramming of DNA replication timing in induced pluripotent stem cells.

Induced pluripotent stem cells (iPSCs) are the foundation of cell therapy. Differences in gene expression, DNA methylation, and chromatin conformation, which could affect differentiation capacity, have been identified between iPSCs and embryonic stem cells (ESCs). Less is known about whether DNA replication timing, a process linked to both genome regulation and genome stability, is efficiently reprogrammed to the embryonic state. To answer this, we compare genome-wide replication timing between ESCs, iPSCs, and cells reprogrammed by somatic cell nuclear transfer (NT-ESCs). While NT-ESCs replicate their DNA in a manner indistinguishable from ESCs, a subset of iPSCs exhibits delayed replication at heterochromatic regions containing genes downregulated in iPSCs with incompletely reprogrammed DNA methylation. DNA replication delays are not the result of gene expression or DNA methylation aberrations and persist after cells differentiate to neuronal precursors. Thus, DNA replication timing can be resistant to reprogramming and influence the quality of iPSCs.

Erratum: The landscape of somatic mutations in lymphoblastoid cell lines.

[This corrects the article DOI: 10.1016/j.xgen.2023.100305.].

Incomplete Reprogramming of DNA Replication Timing in Induced Pluripotent Stem Cells.

Induced pluripotent stem cells (iPSC) are a widely used cell system and a foundation for cell therapy. Differences in gene expression, DNA methylation, and chromatin conformation, which have the potential to affect differentiation capacity, have been identified between iPSCs and embryonic stem cells (ESCs). Less is known about whether DNA replication timing - a process linked to both genome regulation and genome stability - is efficiently reprogrammed to the embryonic state. To answer this, we profiled and compared genome-wide replication timing between ESCs, iPSCs, and cells reprogrammed by somatic cell nuclear transfer (NT-ESCs). While NT-ESCs replicated their DNA in a manner indistinguishable from ESCs, a subset of iPSCs exhibit delayed replication at heterochromatic regions containing genes downregulated in iPSC with incompletely reprogrammed DNA methylation. DNA replication delays were not the result of gene expression and DNA methylation aberrations and persisted after differentiating cells to neuronal precursors. Thus, DNA replication timing can be resistant to reprogramming and lead to undesirable phenotypes in iPSCs, establishing it as an important genomic feature to consider when evaluating iPSC lines.

The landscape of somatic mutations in lymphoblastoid cell lines.

Somatic mutations have important biological ramifications while exerting substantial rate, type, and genomic location heterogeneity. Yet, their sporadic occurrence makes them difficult to study at scale and across individuals. Lymphoblastoid cell lines (LCLs), a model system for human population and functional genomics, harbor large numbers of somatic mutations and have been extensively genotyped. By comparing 1,662 LCLs, we report that the mutational landscape of the genome varies across individuals in terms of the number of mutations, their genomic locations, and their spectra; this variation may itself be modulated by somatic trans-acting mutations. Mutations attributed to the translesion DNA polymerase η follow two different modes of formation, with one mode accounting for the hypermutability of the inactive X chromosome. Nonetheless, the distribution of mutations along the inactive X chromosome appears to follow an epigenetic memory of the active form.

Cell-type specificity of the human mutation landscape with respect to DNA replication dynamics.

The patterns of genomic mutations are associated with various genomic features, most notably late replication timing, yet it remains contested which mutation types and signatures relate to DNA replication dynamics and to what extent. Here, we perform high-resolution comparisons of mutational landscapes between lymphoblastoid cell lines, chronic lymphocytic leukemia tumors, and three colon adenocarcinoma cell lines, including two with mismatch repair deficiency. Using cell-type-matched replication timing profiles, we demonstrate that mutation rates exhibit heterogeneous replication timing associations among cell types. This cell-type heterogeneity extends to the underlying mutational pathways, as mutational signatures show inconsistent replication timing bias between cell types. Moreover, replicative strand asymmetries exhibit similar cell-type specificity, albeit with different relationships to replication timing than mutation rates. Overall, we reveal an underappreciated complexity and cell-type specificity of mutational pathways and their relationship to replication timing.

Generation, genomic characterization, and differentiation of triploid human embryonic stem cells.

Humans are diploid organisms, and triploidy in human embryos is responsible for ∼10% of spontaneous miscarriages. Surprisingly, some pregnancies proceed to triploid newborns that suffer from many neuro-developmental disorders. To investigate the impact of triploidy on human development, we generate triploid human embryonic stem cells (hESCs) by fusing isogenic haploid and diploid hESCs. Comparison of the transcriptome, methylome, and genome-wide replication timing shows general similarity between diploid and triploid hESCs. However, triploid cells have a larger volume than diploid cells, demonstrating decreased surface-area-to-volume ratio. This leads to a significant downregulation of cell surface ion channel genes, which are more essential in neural progenitors than in undifferentiated cells, leading to inhibition of differentiation, and it affects the neuronal differentiation ability of triploid hESCs, both in vitro and in vivo. Notably, our research establishes a platform to study triploidy in humans and points to their pathology as observed in triploid embryos.

The evolution of the human DNA replication timing program.

DNA is replicated according to a defined spatiotemporal program that is linked to both gene regulation and genome stability. The evolutionary forces that have shaped replication timing programs in eukaryotic species are largely unknown. Here, we studied the molecular causes and consequences of replication timing evolution across 94 humans, 95 chimpanzees, and 23 rhesus macaques. Replication timing differences recapitulated the species' phylogenetic tree, suggesting continuous evolution of the DNA replication timing program in primates. Hundreds of genomic regions had significant replication timing variation between humans and chimpanzees, of which 66 showed advances in replication origin firing in humans, while 57 were delayed. Genes overlapping these regions displayed correlated changes in expression levels and chromatin structure. Many human-chimpanzee variants also exhibited interindividual replication timing variation, pointing to ongoing evolution of replication timing at these loci. Association of replication timing variation with genetic variation revealed that DNA sequence evolution can explain replication timing variation between species. Taken together, DNA replication timing shows substantial and ongoing evolution in the human lineage that is driven by sequence alterations and could impact regulatory evolution at specific genomic sites.

DNA replication in early mammalian embryos is patterned, predisposing lamina-associated regions to fragility.

DNA replication in differentiated cells follows a defined program, but when and how it is established during mammalian development is not known. Here we show using single-cell sequencing, that both bovine and mouse cleavage stage embryos progress through S-phase in a defined pattern. Late replicating regions are associated with the nuclear lamina from the first cell cycle after fertilization, and contain few active origins, and few but long genes. Chromosome breaks, which form spontaneously in bovine embryos at sites concordant with human embryos, preferentially locate to late replicating regions. In mice, late replicating regions show enhanced fragility due to a sparsity of dormant origins that can be activated under conditions of replication stress. This pattern predisposes regions with long neuronal genes to fragility and genetic change prior to segregation of soma and germ line. Our studies show that the formation of early and late replicating regions is among the first layers of epigenetic regulation established on the mammalian genome after fertilization.

Genomic signature of Fanconi anaemia DNA repair pathway deficiency in cancer.

Fanconi anaemia (FA), a model syndrome of genome instability, is caused by a deficiency in DNA interstrand crosslink repair resulting in chromosome breakage1-3. The FA repair pathway protects against endogenous and exogenous carcinogenic aldehydes4-7. Individuals with FA are hundreds to thousands fold more likely to develop head and neck (HNSCC), oesophageal and anogenital squamous cell carcinomas8 (SCCs). Molecular studies of SCCs from individuals with FA (FA SCCs) are limited, and it is unclear how FA SCCs relate to sporadic HNSCCs primarily driven by tobacco and alcohol exposure or infection with human papillomavirus9 (HPV). Here, by sequencing genomes and exomes of FA SCCs, we demonstrate that the primary genomic signature of FA repair deficiency is the presence of high numbers of structural variants. Structural variants are enriched for small deletions, unbalanced translocations and fold-back inversions, and are often connected, thereby forming complex rearrangements. They arise in the context of TP53 loss, but not in the context of HPV infection, and lead to somatic copy-number alterations of HNSCC driver genes. We further show that FA pathway deficiency may lead to epithelial-to-mesenchymal transition and enhanced keratinocyte-intrinsic inflammatory signalling, which would contribute to the aggressive nature of FA SCCs. We propose that the genomic instability in sporadic HPV-negative HNSCC may arise as a result of the FA repair pathway being overwhelmed by DNA interstrand crosslink damage caused by alcohol and tobacco-derived aldehydes, making FA SCC a powerful model to study tumorigenesis resulting from DNA-crosslinking damage.

Error-prone repair of stalled replication forks drives mutagenesis and loss of heterozygosity in haploinsufficient BRCA1 cells.

Germline mutations in the BRCA genes are associated with a higher risk of carcinogenesis, which is linked to an increased mutation rate and loss of the second unaffected BRCA allele (loss of heterozygosity, LOH). However, the mechanisms triggering mutagenesis are not clearly understood. The BRCA genes contain high numbers of repetitive DNA sequences. We detected replication forks stalling, DNA breaks, and deletions at these sites in haploinsufficient BRCA cells, thus identifying the BRCA genes as fragile sites. Next, we found that stalled forks are repaired by error-prone pathways, such as microhomology-mediated break-induced replication (MMBIR) in haploinsufficient BRCA1 breast epithelial cells. We detected MMBIR mutations in BRCA1 tumor cells and noticed deletions-insertions (>50 bp) at the BRCA1 genes in BRCA1 patients. Altogether, these results suggest that under stress, error-prone repair of stalled forks is upregulated and induces mutations, including complex genomic rearrangements at the BRCA genes (LOH), in haploinsufficient BRCA1 cells.

Replication stress impairs chromosome segregation and preimplantation development in human embryos.

Human cleavage-stage embryos frequently acquire chromosomal aneuploidies during mitosis due to unknown mechanisms. Here, we show that S phase at the 1-cell stage shows replication fork stalling, low fork speed, and DNA synthesis extending into G2 phase. DNA damage foci consistent with collapsed replication forks, DSBs, and incomplete replication form in G2 in an ATR- and MRE11-dependent manner, followed by spontaneous chromosome breakage and segmental aneuploidies. Entry into mitosis with incomplete replication results in chromosome breakage, whole and segmental chromosome errors, micronucleation, chromosome fragmentation, and poor embryo quality. Sites of spontaneous chromosome breakage are concordant with sites of DNA synthesis in G2 phase, locating to gene-poor regions with long neural genes, which are transcriptionally silent at this stage of development. Thus, DNA replication stress in mammalian preimplantation embryos predisposes gene-poor regions to fragility, and in particular in the human embryo, to the formation of aneuploidies, impairing developmental potential.

Telomere-to-telomere human DNA replication timing profiles.

The spatiotemporal organization of DNA replication produces a highly robust and reproducible replication timing profile. Sequencing-based methods for assaying replication timing genome-wide have become commonplace, but regions of high repeat content in the human genome have remained refractory to analysis. Here, we report the first nearly-gapless telomere-to-telomere replication timing profiles in human, using the T2T-CHM13 genome assembly and sequencing data for five cell lines. We find that replication timing can be successfully assayed in centromeres and large blocks of heterochromatin. Centromeric regions replicate in mid-to-late S-phase and contain replication-timing peaks at a similar density to other genomic regions, while distinct families of heterochromatic satellite DNA differ in their bias for replicating in late S-phase. The high degree of consistency in centromeric replication timing across chromosomes within each cell line prompts further investigation into the mechanisms dictating that some cell lines replicate their centromeres earlier than others, and what the consequences of this variation are.

High-throughput analysis of single human cells reveals the complex nature of DNA replication timing control.

DNA replication initiates from replication origins firing throughout S phase. Debate remains about whether origins are a fixed set of loci, or a loose agglomeration of potential sites used stochastically in individual cells, and about how consistent their firing time is. We develop an approach to profile DNA replication from whole-genome sequencing of thousands of single cells, which includes in silico flow cytometry, a method for discriminating replicating and non-replicating cells. Using two microfluidic platforms, we analyze up to 2437 replicating cells from a single sample. The resolution and scale of the data allow focused analysis of replication initiation sites, demonstrating that most occur in confined genomic regions. While initiation order is remarkably similar across cells, we unexpectedly identify several subtypes of initiation regions in late-replicating regions. Taken together, high throughput, high resolution sequencing of individual cells reveals previously underappreciated variability in replication initiation and progression.

Comprehensive analysis of DNA replication timing across 184 cell lines suggests a role for MCM10 in replication timing regulation.

Cellular proliferation depends on the accurate and timely replication of the genome. Several genetic diseases are caused by mutations in key DNA replication genes; however, it remains unclear whether these genes influence the normal program of DNA replication timing. Similarly, the factors that regulate DNA replication dynamics are poorly understood. To systematically identify trans-acting modulators of replication timing, we profiled replication in 184 cell lines from three cell types, encompassing 60 different gene knockouts or genetic diseases. Through a rigorous approach that considers the background variability of replication timing, we concluded that most samples displayed normal replication timing. However, mutations in two genes showed consistently abnormal replication timing. The first gene was RIF1, a known modulator of replication timing. The second was MCM10, a highly conserved member of the pre-replication complex. Cells from a single patient carrying MCM10 mutations demonstrated replication timing variability comprising 46% of the genome and at different locations than RIF1 knockouts. Replication timing alterations in the mutated MCM10 cells were predominantly comprised of replication delays and initiation site gains and losses. Taken together, this study demonstrates the remarkable robustness of the human replication timing program and reveals MCM10 as a novel candidate modulator of DNA replication timing.

Integrated genome and transcriptome analyses reveal the mechanism of genome instability in ataxia with oculomotor apraxia 2.

Mutations in the SETX gene, which encodes Senataxin, are associated with the progressive neurodegenerative diseases ataxia with oculomotor apraxia 2 (AOA2) and amyotrophic lateral sclerosis 4 (ALS4). To identify the causal defect in AOA2, patient-derived cells and SETX knockouts (human and mouse) were analyzed using integrated genomic and transcriptomic approaches. A genome-wide increase in chromosome instability (gains and losses) within genes and at chromosome fragile sites was observed, resulting in changes to gene-expression profiles. Transcription stress near promoters correlated with high GCskew and the accumulation of R-loops at promoter-proximal regions, which localized with chromosomal regions where gains and losses were observed. In the absence of Senataxin, the Cockayne syndrome protein CSB was required for the recruitment of the transcription-coupled repair endonucleases (XPG and XPF) and RAD52 recombination protein to target and resolve transcription bubbles containing R-loops, leading to genomic instability. These results show that transcription stress is an important contributor to SETX mutation-associated chromosome fragility and AOA2.

Delayed DNA replication in haploid human embryonic stem cells.

Haploid human embryonic stem cells (ESCs) provide a powerful genetic system but diploidize at high rates. We hypothesized that diploidization results from aberrant DNA replication. To test this, we profiled DNA replication timing in isogenic haploid and diploid ESCs. The greatest difference was the earlier replication of the X Chromosome in haploids, consistent with the lack of X-Chromosome inactivation. We also identified 21 autosomal regions that had delayed replication in haploids, extending beyond the normal S phase and into G2/M. Haploid-delays comprised a unique set of quiescent genomic regions that are also underreplicated in polyploid placental cells. The same delays were observed in female ESCs with two active X Chromosomes, suggesting that increased X-Chromosome dosage may cause delayed autosomal replication. We propose that incomplete replication at the onset of mitosis could prevent cell division and result in re-entry into the cell cycle and whole genome duplication.

The genetic architecture of DNA replication timing in human pluripotent stem cells.

DNA replication follows a strict spatiotemporal program that intersects with chromatin structure but has a poorly understood genetic basis. To systematically identify genetic regulators of replication timing, we exploited inter-individual variation in human pluripotent stem cells from 349 individuals. We show that the human genome's replication program is broadly encoded in DNA and identify 1,617 cis-acting replication timing quantitative trait loci (rtQTLs) - sequence determinants of replication initiation. rtQTLs function individually, or in combinations of proximal and distal regulators, and are enriched at sites of histone H3 trimethylation of lysines 4, 9, and 36 together with histone hyperacetylation. H3 trimethylation marks are individually repressive yet synergistically associate with early replication. We identify pluripotency-related transcription factors and boundary elements as positive and negative regulators of replication timing, respectively. Taken together, human replication timing is controlled by a multi-layered mechanism with dozens of effectors working combinatorially and following principles analogous to transcription regulation.

Replication timing analysis in polyploid cells reveals Rif1 uses multiple mechanisms to promote underreplication in Drosophila.

Regulation of DNA replication and copy number is necessary to promote genome stability and maintain cell and tissue function. DNA replication is regulated temporally in a process known as replication timing (RT). Rap1-interacting factor 1 (Rif1) is a key regulator of RT and has a critical function in copy number control in polyploid cells. Previously, we demonstrated that Rif1 functions with SUUR to inhibit replication fork progression and promote underreplication (UR) of specific genomic regions. How Rif1-dependent control of RT factors into its ability to promote UR is unknown. By applying a computational approach to measure RT in Drosophila polyploid cells, we show that SUUR and Rif1 have differential roles in controlling UR and RT. Our findings reveal that Rif1 acts to promote late replication, which is necessary for SUUR-dependent underreplication. Our work provides new insight into the process of UR and its links to RT.

Cancer Mutational Processes Vary in Their Association with Replication Timing and Chromatin Accessibility.

Cancer somatic mutations are the product of multiple mutational and repair processes, both of which are tightly associated with DNA replication. Distinctive patterns of somatic mutation accumulation, termed mutational signatures, are indicative of processes sustained within tumors. However, the association of various mutational processes with replication timing (RT) remains an open question. In this study, we systematically analyzed the mutational landscape of 2,787 tumors from 32 tumor types separately for early and late replicating regions using sequence context normalization and chromatin data to account for sequence and chromatin accessibility differences. To account for sequence differences between various genomic regions, an artificial genome-based approach was developed to expand the signature analyses to doublet base substitutions and small insertions and deletions. The association of mutational processes and RT was signature specific: Some signatures were associated with early or late replication (such as SBS7b and SBS7a, respectively), and others had no association. Most associations existed even after normalizing for genome accessibility. A focused mutational signature identification approach was also developed that uses RT information to improve signature identification; this approach found that SBS16, which is biased toward early replication, is strongly associated with better survival rates in liver cancer. Overall, this novel and comprehensive approach provides a better understanding of the etiology of mutational signatures, which may lead to improved cancer prevention, diagnosis, and treatment. SIGNIFICANCE: Many mutational processes associate with early or late replication timing regions independently of chromatin accessibility, enabling development of a focused identification approach to improve mutational signature detection.